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Abstract 

We consider a setup where R-parity is violated in the framework of spht 
supersymmetry. The out-of-equihbrium decays of heavy squarks successfuhy 
lead to the generation of a baryon asymmetry. We restrict the R-parity vio- 
lating couplings to the baryon number violating subset to keep the neutralino 
sufficiently stable to provide the dark matter. The observed baryon asymme- 
try can be generated for squark masses larger than 10^^ GeV, while neutralino 
dark matter induces a stronger bound of 10^^ GeV. Some mass splitting be- 
tween left- and right-handed squarks may be needed to satisfy also constraints 
from gluino cosmology. 
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1 Introduction 



The general supersymmetric extension of the standard model (SM) contains baryon 
and lepton number violating interactions [1,2]. Often a discrete symmetry, called 
R-parity, is imposed to eliminate these operators. In conventional supersymmet- 
ric models there are stringent bounds on many of the R-parity violating couplings 
induced by proton decay, neutrino masses, flavor violation, etc. The breaking of R- 
parity also has important cosmological implications. The baryon and lepton number 
violating interactions may erase a baryon asymmetry created in the early universe. 
They also induce the decay of the hghtest supersymmetric particle (LSP). Only in 
the case of extremely small R-parity violating couplings does the LSP remain a 
viable dark matter candidate. 

Most of these constraints are considerably relaxed in the context of split super- 
symmetry, where all scalars, except for a single (and finely tuned) Higgs boson, have 
masses much larger than the electroweak scale, rh < 10^^ GeV [3-5]. For instance, 
bounds from proton decay have been discussed in ref. [6]. Taking rh — 10^^ GeV, 
the product of baryon and lepton number violating couplings can be as large as 
10~^. In conventional supersymmetry it is bounded to be below 10~^^ to guarantee 
a sufficiently long-lived proton. 

The lepton number violating superpotential contains bilinear terms, /i^LjiJ2- 
They lead to a tree-level mixing between the neutrinos and neutralinos which is 
not suppressed by the heavy mass scale of split supersymmetry. This mixing gen- 
erates a neutrino mass at the tree-level which constrains the /i'-terms to be smaller 
than about 1 MeV [7]. It also induces the decay of the lightest neutralino, e.g. into 
a neutrino and a (virtual) Z boson. Keeping the neutralino lifetime larger than 
the age of the universe requires < 10~^ eV. In contrast, neutralino decay via R- 
parity violating Yukawa couplings (trilinear terms) is indeed suppressed by the large 
sfermion masses [8] . The lepton number violating Yukawa couplings induce //'-terms 
in the low energy effective action by quantum corrections even if they are set to zero 
at the tree-level [9] . Thus, an enormous amount of tuning is needed in this setting 
to keep the lightest neutrahno sufficiently stable to provide the dark matter. 

In this article we show how the baryon asymmetry of the universe could be 
generated by the out-of-equilibrium decays of heavy sfermions in the framework 
of split supersymmetry^. In this context, neutralino dark matter is an important 
hnk to tie the masses of the gauginos and higgsinos to the electroweak scale [3-5] . 
To avoid rapid neutralino decay by lepton number violating bilinears, we restrict 
R-parity violation to the baryon number violating operators. 

There have been attempts to use R-parity violating interactions to generate the 
baryon asymmetry in weak scale supersymmetry [11]. However, these typically re- 
quire some additional ingredients to provide the necessary departure from thermal 
equilibrium, such as non-thermal production of squarks at the end of inflation, late 
decay of gra\itiiios or axinos. etc. We will show that in split supersymmetry ther- 

^See ref. [10] for an alternative proposal. 
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mal baryogenesis is successful for squark masses larger than about 10^^ GeV. As 
it turns out, the ncutralino lifetime induces a somewhat stronger constraint on the 
squark masses of about 10^^ GeV. Gluino cosmology may require some mass splitting 
between left- and right-handed squarks. 

The organization of the paper is as follows. In sec. 2 we discuss how a CP- 
asymmetry arises in the squark decays. It requires the interference of tree-level 
and two-loop diagrams. In sec. 3 we write down the Boltzmann equations, which 
describe the baryon production process. Their numerical solution together with 
some analytical approximations will be presented in sec. 4. Constraints induced by 
the gluino and neutralino hfetimes will be discussed in sec. 5. In sec. 6 we present 
our conclusions. 

2 CP-asymmetry in squark decays 

According to the discussion in the introduction we consider the following R-parity 
violating interactions in the superpotential [1,2]^ 

W^, = KjkU^Dpi, (1) 

where Xijk = —\kj- These operators violate baryon number but not lepton number. 
The proton is therefore stable. The form of the superpotential (1) could, for instance, 
be guaranteed by lepton parity. 

Since A-terms, gaugino masses and the /x-parameter are small compared to the 
split supersymmetry scale rh [3-5], scalar masses for Uf and are the only soft 
terms relevant to our discussion. We can diagonalize them by supersymmetric ro- 
tations. Then flavor transitions can only be induced by the R-parity violating cou- 
plings (1) or by ordinary Yukawa couplings. We make the (conservative) assumption 
that the effects of the latter are small, and ignore them. 

We assume a (somewhat) hierarchical spectrum of the SU(2) singlet ("right- 
handed") squarks"^. Then baryogenesis is dominated by the decay of the lightest of 
these states, which we take to be the right-handed up squark u. It will turn out that 
this choice simplifies the Boltzmann equations, which describe the baryogenesis pro- 
cess. Taking another right-handed squark would change the final baryon asymmetry 
only by a factor of order 1. 

The spectrum of the SU(2) doublet ("left-handed") squarks is not directly re- 
lated to the baryogenesis process. If some of them are lighter than u, they only 
enter the Boltzmann equations as additional degrees of freedom, which carry part 
of the baryon number. Later on this possibility will turn out to be helpful to rec- 
oncile the constraints induced by the gluino and neutralino lifetimes. A somewhat 

^Note that in the Uterature these couphngs arc usuaUy denoted by A^^j,. 

"^For the SU(2) singlet states we use the notation u = etc., so that they are counted with 

positive baryon number. 
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lighter left-handed squark can speed up the gluino decay through supergauge in- 
teractions, making it more easy to satisfy constraints from gluino cosmology [12]. 
The neutralino decay by R-parity violating couplings is still governed by the (heav- 
ier) right-handed squarks. Thus the neutralino can remain sufficiently long-lived to 
provide the dark matter of the universe. 

At least four of the nine couplings of the superpotential (1) have to be present 
to allow for CP violation in the up squark decays. In the following wc assume two 
of them to be A112 and A123. Two other couplings have to be related to a different, 
heavier squark, which we take to be t. The first two operators induce baryon number 
violating decays of the up squark. At tree-level the partial decay widths are 

r{u^ d2dz) = 823^^0 

where To denotes the total decay width and mn the up squark mass. Note that the 
d quarks are also right-handed (SU(2) singlet) states. For simplicity we ignore a 
possible decay into the dids channel. The up squark also decays by the supergauge 
interaction at a rate 

2 

r{u ug) = EgVo = -^asnnu. (3) 

Finally, the up squark can decay through its Yukawa coupling This contribution 
is suppressed by the small value of Uu and ignored in the following. 

At the loop-level, the decay widths for u and and its antiparticle u* become 
different in the presence of CP- violation 

v{u^d^d2) = rc(si2-^) 

T{u d2d3) = Yd (^^23 - 
r(u^ug) = r^(i?, + |). (4) 

Going to antiparticles, the sign of the CP- violating contributions are reversed. Since 
u and its anti-particle u* have the same width, 

eg = ei2 + £23 (5) 

holds. 

A CP-asymmetry in the decay width is generated by the interference of tree-level 
and 2-loop amplitudes^. The diagram of fig. la leads to 

ei2 oc lm(Aii2Ai23A3i2A323) oc sm{6), (6) 
^Note that 1-loop contributions are suppressed by the smaU gaugino masses and A-parameters. 
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Figure 1: 2- loop contribution (a) to €±2 and its approximation by a 1-loop diagram 
(b). 

where A312 and A323 arc related to the (heavier) top squark, and 6 defines the effective 
CP-phase. Instead of computing the 2-loop diagram, we approximate the effect of 
the stop by a non-diagonal gluino vertex 

y/2e^zgfgdld^. (7) 

Computing the remaining 1-loop diagram of fig. lb , we obtain for the CP-asymmetry 

2 

ei2 = -Q;sSin(5)Si26'i3 

The loop function is given by 

/(x) = 2^1-iln(l + a;)j aix+C>(x2). (9) 

Of course, there are more diagrams that contribute^, so this only a rough estimate. 
Note that 612 can be enhanced by taking the ratio A123/A112 to be large. The CP- 
asymmetry 623 in the d2dz channel is obtained hy 5 ^ —5 and — > mj . 
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3 Boltzmann equations 

The evolution of the baryon asymmetry is governed by a set of Boltzmann equations. 
We assume that kinetic equilibrium is maintained and approximate the phase-space 
densities by Maxwell-Boltzmann distributions, 

f/E) = e-^^-^^y^, (10) 

^The off-diagonal gluino vertex can also be generated by ordinary Yukawa couplings, with left- 
handed squarks running in the loop. 
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where /li is the chemical potential of the ith particle species. It describes the devi- 
ation from chemical equilibrium. 

3.1 Up squark decay 

The number density of right-handed up squarks is governed by the Boltzmann equa- 
tion 

dn. 



dt 



+ 2>Hn^ = -Da - Di2 - D 



'23- 



11) 



Here H = 1.66^^7^^ denotes the Hubble parameter, with Mp = 1.22 x lO^'^ 
GeV. If all scalars are heavier than u, the number of degrees of freedom in the plasma 
is = 131.25. For every left-handed squark doublet, which is light, increases by 
12. To first order in the CP-asymmetries, the collision terms for the different decay 
channels are given by 



Du = 
D23 = 



7d 
Id 



(i?i2 - f ) - nf{l - U - U) {b,2 + 'f) 



(12) 



where C,i = fJ-i/T. To obtain this result wc approximated cxp(,^j) ^ 1 + and used 
= — for the antiquark chemical potentials. The quark chemical potentials are 
first order in the CP-asymmetries, so that to first order we can drop terms like ^e. 
We used energy conservation to eliminate the quark distribution functions and CPT 
invariance to obtain the rates for the inverse decay processes. The equilibrium num- 
ber density n'j^ is computed from the distribution function with vanishing chemical 
potential and includes three degrees of freedom. The thermal decay rate averaged 
with equilibrium distribution function is [13] 



where Ki and K2 are Bessel functions and 

rriu 



z — 



T 



(13) 



(14) 



The evolution equation for u* is obtained by reversing the signs of the CP- violating 
contributions and replacing the quark chemical potentials by that of antiquarks, and 

vice versa. 

It is useful to write the number densities of u and u* as [13] 



(15) 
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Up to linear order in the CP-violating asymmetries, we obtain 



^ + 3HN^ = -^n{N^ - K^) (16) 

The equilibrium density = 2rf^ = {3T^/tt'^)z^K2{z) contains six degrees of 
freedom. 

We normahze the number densities to the entropy density s — (27r^5i* 7^/45), 
where is again the number of degrees of freedom in the plasma, such that 

TV- 

Y, = ^. (17) 
s 

Normalization by s removes the term proportional to the Hubble parameter from 
the Boltzmann equation. It is convenient to transform to the dimensionless variable 
z, using dY/dt — {z/H{T — mu)){dY / dz), to arrive at 



dY+ 
dz 

where 



-l}K(Y^-Y^), (18) 



Assuming that the decay rate is dominated by the gluino channel, we have K k, l.bx 
10^ (10^^ GeV/rriu), where the running strong coupling is as{10^^ GeV) = 0.035. A 
large value of K signals that the up squark decay is close to equilibrium. Along the 
same lines we arrive at 

dY 

= -PK{Y_ - - B,2{U + U) - B23{U + U)])- (20) 

We see that the inverse decays lead to a back reaction of the quark densities on the 
evolution of Y_. 



3.2 Evolution equations for the quark densities 

Before computing the evolution of the quark densities we have to check what other 
than the R-parity violating interactions are relevant at temperatures around lO^^GeV 
which we are going to consider. 

i) Strong sphalcrons transform left-handed quarks into right-handed ones. Using the 
results of ref. [14] and taking into account the running of as, we find these processes 
to be in equilibrium for T < 5 x 10^^ GeV (see also ref. [15]). 

ii) Weak sphalerons [16] become fast at temperatures T <2 x 10^^ GeV. 
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iii) The scattering rate of Yukawa interactions has been estimated as Ty ^ 0.02y^T 
[17]. The top Yukawa couphng therefore is in cquiUbrium for T < 10^^ GeV. At 
T < 1 X 10^^ GeV also the bottom and tau Yukawa interactions become relevant. 

In the following we will assume the strong sphalerons and the top Yukawa inter- 
action to be in equilibrium. All other Yukawa couplings and the weak sphalerons will 
be neglected. This approximation is vahd in the temperature range 1 x 10^^ GeV < T < 5 x 
10^^ GeV and leaves the leptons entirely out of the baryon generation process. 

The relevant operators are 

gu*u, udid2, ud2d3, gd^di, 

uctdxd2d^qiqiq2q2M'i- (21) 

The quark species c, qi and q2 are only produced by strong sphalerons. Therefore 
the number densities are related as 



2N, = -N,, = -Nq, = 2Ng. (22) 

Here we use the number densities of particles minus antiparticles, i.e. Nc = ric — 
ncCtc. The operators (21) conserve four independent U(l) charges, which imply the 
constraints 

Nt + Ng,+Ng = 

Nd, - Nd, + Nd, - N, = 
N_ + Nu-Nd, = 

Nt + Nn-Ng = 0. (23) 

Taking the top Yukawa and strong sphaleron interactions to be in equilibrium forces 
the corresponding linear combinations of chemical potentials to vanish. What these 
relations mean in terms of number densities depends on the number of degrees 
of freedom (d.o.f.) related to each chemical potential. To be specific, we assume 
that one Higgs doublet is much heavier than rriu and the supergauge interactions 
maintain the equilibrium between the Higgs and higgsino densities. The Higgs 
chemical potential then represents 12 d.o.f., like q^. As already mentioned, we 
can allow for some left-handed squarks to be light, which we take to be qi. Again 
supergauge interactions will equilibrate the and qi densities and effectively double 
the d.o.f. of qi. The equilibrium conditions then take the form 

-2Nt + Ng, + NH = 
N^ + Nt + Na,+Na, + Nd,-N,,+ (^3 + j^Ng = 0, (24) 

where k — 1 (2) if qi is heavy (light). Using eqs. (22), (23) and (24) we can express 
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all quantities in terms of N^, Nu and N^^ as 

Nd, = N_ + Nu 



2k 3k 
" " 5k + 2 5A; + 2 " 

Nt = 0. (25) 

We can use these relations to express baryon number as 

Nb = ^{N_ + Nu + N, + Nt + Nd, + Nd,+Nd, + Ng, + Ng,+Ng,) 

= Z^7V_ + ^7V„. (26) 

Note that A^^^ does not enter here directly. This is an accident of the approximations 
we arc using. For instance, it would no longer hold if the bottom Yukawa coupling 
were taken into account. 

To compute the generated baryon asymmetry we need two evolution equations 
in addition to eqs. (18), (20). We take them to be 

^(y„-n) = Dg-S^2-S2s-T,2-T2s-S[^-S!,^-Ti^-Ti^-CP conj. 

^(Fd, -Fe) = -Du-Su-Tu-S,s-S[^-T,s-T[^-CF coni. (27) 

Subtracting Yc in these equations removes the strong sphaleron rate from the right- 
hand side, which we have taken to be in equilibrium. The inverse decays D12, D23 
and Dg are given by eq. (12). We can use Uu — riu ~ 2^un^'^, etc., with rf^ = ST^/tt^, 
to replace the chemical potentials by particle densities, and we define 

A{z) = -^ = -K,{z). (28) 
rig z 

The scatterings induced by u exchange are 

-A{z)BgBu{Y^ + Yd, + Yd,) 

S23-S23 = {Yu + Yd,+Yd,hs,s+Y;''^n{Bge23 + B2seg) 
-A{z)BgB2s{Yu + Yd,+Yd,) 

Si3 - S^3 = (Yd, - Yd^y/s,, + Y;^Jn{B23e^2 - ^12623) - A{z)B^2B23{Yd, - Yd,) 

T12-T12 = 2{Y^ + Yd,+Yd,)-fT,, 

T23-T23 = 2{Y^ + Yd,+Yd,hT,s 

Ti3-ri3 = 2{Yd,-Yd,)^T,s- (29) 
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Here 5'i2, S23 and S13 correspond to the processes tg <-> ^1^2, tg <-> c?2C^3 and 
did2 ^ a?2C^3- Note that for these s-channel processes the resonant parts have been 
subtracted since they are aheady included in the Bohzmann equations [13]. T12, T23 
and Ti3 denote tdi^2 ^ 5<^2,i, ^<^2,3 ^ 5<^3,2 and did2,3 <i3,2<^2- Also the exchange 
of ^3 contributes to the above mentioned scatterings 

S23-s'23 = (y„ + y,, + y,3)75^3 
T{2-T{2 = 2(i; + r,,+y,j7T^3 

7^23-^23 = 2(F„ + F,, + r,3)7T^3. (30) 

In an approximation we add these contributions incoherently, which is valid when 
either u or ^3 exchange dominates. There are also contributions to the scatterings 
by exchange of t, which we neglect. 

For temperatures considerably smaller than the masses of the exchanged squarks 
we obtain for the scatterings 

352 2 1 



rUu 352 2 1 
ISr + 27Tfi = ^3^^i2"^^(-^x + >di + Yd^) 



75i3 + 27ti3 = ^16^ 112-^113^^ ^1 ~ 

^S', + 2^T', = '^^jfsXl,,a.\^{Y^ + Y,, + Y,,). (31) 
H 37r^ m*~ 
03 

Note that 75^ + 27^^^ is suppressed by a large factor m\/mj-^ ~ 1/100 compared to 
'^Sa + 27ta ^'^d can be neglected in the following. 

Now we can use the relations (25) to eliminate Yd^ and Y^,^ in the inverse decay 
and scattering rates to turn eqs. (20), (27) into a closed set of equations for y_, Yu 
and Yd^ . The result is given in eq. (44) in the appendix. 



4 Analytical and numerical solutions 

To understand some features of the numerical solution of the Boltzmann equations 
(44), let us first discuss some analytical approximations. 

The departure from equihbrium in the evolution of 1+ is of order 1/K. Since for 
rriu < 10^^ GeV we always have X » 1, we can expand in powers of 1/K (see, for 
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instance ref. [18]) To leading order we obtain 



1 dYl'' 
~PK dz 

— y^i 
Kz + 

135 z 135 ^/z 

Non-equilibrium in the CP- violating densities is governed by some effective value of 
which involves the baryon number violating couplings. 

Let us start with the somewhat simpler case Byi = B23. Then the down quark 
number densities combine as Frf^+2Fd2+^d3) since ■Jsi2,Ti2 = 7523,T23- (This is not true 
for the scatterings induced by d^, which however are suppressed by m\/rnh ~ 1/100 
and can be neglected.) As a result, we can form a closed set of equations for y_ and 
Yb-i which is given in the appendix (45). In the limit of large the equation for 
y_ reduces to 

We also have neglected terms higher order in A[z), which are exponentially small 
at late times, i.e. for 2; ^ 1. We end up with a single equation for the baryon 
asymmetry 



dYs 
dz 



= (5k\^-^{Y^-Y^)-A{z)B,,{1-2B,,)^-^±^Yb 
llA; + 5 m^i352 2 1 



4)t + l ^i/ 37r2^^ii2"^z4- 

We observe that the effective value of K is reduced proportionally to B12 and depends 
somewhat on k. 

For large values of B12K we can now solve for Yb using the method of steepest 
descent. Assuming that the inverse decays dominate, the freeze-out value of z follows 
from 



5/2 23/2 4^ + I 



■B,2{1-2B,2)K, (35) 



which leads to a logarithmic dependence of Zf on K. The final baryon asymmetry 
is then given by 



^286^ 135 4fc + l / 4zf 

79 2 7ry^g,ZfB^2{l - 2B^2)K Ilk + 5\J 2zf - 5' ^ ^ 

where we added the weak sphaleron factor of 28/79. Using eq. (8) we observe that 
the leading dependence of rjB on B12 cancels. This means that smaller values of the 
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Figure 2: (a) The baryon asymmetry in units of 10 as a function of A112 for 
ms = 5 X 10^3 QgY^ Bu = B23, k = 1, 9^ = 0.1 and m^- /m? = 10. The full 
(dashed) line includes (neglects) scatterings. The horizontal dashed lines indicate 
the observed range, (b) The baryon asymmetry with (without) a light qi in solid 
(dashed) line. Scatterings are neglected. 



baryon number violating couplings reduce by the same amount the CP-asymmetry 
and the washout by inverse decays. This compensation works until B12K ~ 1. 
Then the washout is no longer important and the baryon asymmetry goes to zero 
proportional to B12. For very large values of B12K the washout is dominated by 
scatterings. Then zj (x K^^^ and the baryon asymmetry is exponentially damped, 

In the general case of B12 7^ eq. (34) generalizes to 

^ = |^(r+ - Yl^) - '^-A{z)B,2{^ - 7B,2 - 8523)1^^} , (37) 

where we have neglected scatterings and restricted ourselves to €23 = and k = 1. 
Note that the washout term is proportional to Bu, i.e. the effective value of K is 
~ B12K. It can be made small while keeping a large value of B23. This allows 
us to preserve the di part of the baryon asymmetry in the presence of a large A123, 
which can be used to increase the CP asymmetry. A larger value of B23 decreases 
the washout. If €23 is non-zero, its contribution to the baryon asymmetry is washed 
out with f^eff ~ B23K, i.e. it is not protected by a small value B12. 

In the numerical evaluations of eq. (44) we take a maximal CP-phase sin 5 = 1, 
^13 = 0.1 and mj /m| = 10. In fig. 2a we show the baryon asymmetry for B12 = B23 
as a function of A112. The other parameters are k = 1 and mn = 5 x 10^^ GeV. 
The dashed horizontal lines indicate the observed value of the baryon asymmetry, 
riB = (0.89 ± 0.04) x 10"^° [19]. The dashed curve takes into account only the 
inverse decays in the collision terms of eq. (44). In agreement with eq. (36), the 
baryon asymmetry stays almost constant for A112 > 0.02. For smaller values the 
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Figure 3: (a) The baryon asymmetry in units of 10 as function of A112 for rriu = 
1 X 10^3 QgY^ = 0.25, kg^ = 1, = 0.1 and /ml = 10. The full (dashed) 
line includes (neglects) scatterings, (b) The baryon asymmetry with (without) a 
light qi in solid (dashed) line. Scatterings are neglected. 



washout becomes negligible and we have ?7_b oc Bu oc A^^2- The solid line takes into 
account also the scatterings that lead to a strong damping of the baryon asymmetry 
for A112 > 0.2. In fig. 2b we investigate the influence of a light left-handed squark 
qi on the baryon asymmetry, again neglecting scatterings. The difference is quite 
small since two effects partially cancel: the hght gi {k = 2) reduces r]B by adding 
d.o.f. to the plasma. It also reduces the washout, as we observe from eq. (34). For 
smaller up squark masses the maximal baryon asymmetry decreases approximately 
as TjB oc niu. So up squark masses above a few times 10^^ GeV are needed to generate 
the observed baryon asymmetry in the case B12 = B23. 

A larger baryon asymmetry can be produced if A123 > A112. In this case eu is 
enhanced while the washout remains small due to A112. In fig. 3a we take A123 = 0.25, 
rriu = 1 X 10^^ GeV and k = 1. Solid (dashed) lines indicate again that scatterings 
are included (neglected). Now the baryon asymmetry rises proportionally to I/A112. 
For A112 < 0.01 the washout becomes small and tjb oc A112. Fig. 3b shows that the 
inclusion of a light qi does hardly make any difference. In fig. 4 we compare the 
baryon asymmetry for A123 = 0.25 and 0.1. All other parameters are taken as in 
fig. 3a. We observe that rjB scales approximately as A123. 

The value of A112, where the baryon asymmetry becomes maximal, i.e. where the 

— 1/2 

washout becomes ineffective, scales as m- . Using this optimal value of A112, the 
baryon asymmetry is therefore roughly given by 



riB ^ 10"%in((5)A 



The observed baryon asymmetry can therefore be generated for rriu ^ 10^^ GeV.^ 
'''For rriu ^ 10^^ GeV the tau and bottom Yukawa interactions as well as the weak sphalerons 
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Figure 4: The baryon asymmetry in units of 10 as function of A112 for = 
1 X 10^^ GeV, kq^ = 1, 6'i3 = 0.1 and m^- /m\ = 10. Solid (dashed) hues indicate 
A123 = 0.25 (0.1). Scatterings are included. 



For rriu = 1 x 10^'^ GeV the coefficient of the off-diagonal gluino vertex can be chosen 
to be smaller by a factor of 10 than used in the figures. 



5 Long-lived neutralinos and gluinos 

By exchange of the lightest right-handed squark, u, the lightest neutralino predom- 
inantly decays as 

X? ^ Mrf24 (39) 

and the corresponding antiparticles. Again we assumed that i?23 ^ B12, so that the 
decay into udid2 is suppressed. The lifetime associated with this process is 

1 X lO^'^'^yr / Tfi- \ / \^ 
y2^?23 VSOgJ^J Vl0i3GevJ ■ 

where y denotes the coupling 'yc^uu. While for a bino-like Xi one has y = ^\/2gi, 
this coupling is given by the up Yukawa coupling for a higgsino-like LSP, and even 
vanishes for a pure wino LSP. Even with yXi23 ~ 1 the LSP lifetime is expected 
to be much larger than the age of the universe. It can be made even longer lived 
by suppressing y. If we assume gaugino mass unification, the lightest neutralino is 
predominantly a bino in the range of small M2 [4, 20]. Then we estimate y"^ ~ 1/20. 

Stringent constraints on the lifetime of an instable dark matter particle follow 
from the production of antiprotons in its decay (39) [21] 

^xl°">-'(5^)"' («) 

can no longer be neglected. We expect the corresponding change in tje to be at most a factor of 
order unity. 
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Similar constraints can also be derived from positron production [22]. This implies 
that 

It means that for a bino-hke LSP we need rriu > 10^^ GeV to be in agreement 
with the anti-proton bound. Light neutralino masses are clearly favored. In split 
supersymmetry there is the possibility of having a wino- or higgsino as the LSP, 
with a mass of 2.0-2.5 and 1.0-1.2 TeV, respectively [4,20]. In this case y can be 
very much suppressed, but it requires some tuning to keep the bino content at the 
level of 10~^ to compensate for the larger neutrahno masses. It seems to be very 
difficult to use the suppression of y to allow for < 10^^ GeV. Moreover, for very 
small values of y, loop corrections and the contributions from heavier right-handed 
squarks become important. 

A striking signal of split supersymmetry is the very long-lived gluino. Its lifetime 
by R-parity conserving decays has been estimated as [23] 

Heavy isotope searches induce an upper bound on the sfermion mass scale of m < 10^" 
GeV [3]. This is consistent with the lower bound wc previously derived from the 
neutralino lifetime. Taking m(= rriu) = 10^^ GeV and nig = 1 TcV leads to a gluino 
lifetime of about 10^ yr. If B23 ~ 1 and u is the lightest squark, the gluino could 
have a sizable baryon number violating branching ratio. 

Gonstraints from gluino cosmology depend crucially on the gluino annihilation 
cross section after the QGD phase transition. In ref. [12] it has been argued that the 
latter is set by the de Broglie wavelength of the gluino rather than the geometric 
cross section. In this case a stronger bound of m(= m^) = 10^^ GeV can be derived 
from diffuse gamma rays. For nig > 500 GeV an even more restrictive bound is 
induced by big bang nucleosynthesis. To meet these constraints we would have to 
assume a mass splitting between the left- and right-handed squarks. The right- 
handed squarks must be heavier than about 10^^ GeV to keep the neutralino stable 
enough. Some left-handed squarks, e.g. qi have to be lighter than about 10^^ GeV to 
speed up gluino decay. In this case the baryon number violating branching ratio of 
the gluino is highly suppressed by {rhq^/rhuY ■ Gauge couplings unification is hardly 
affected by such a mass splitting. Taking m^^ = 10^^ GeV, the low energy value 
of cts is increased by about 2 x 10~^ with respect to the case where all squarks are 
degenerate at m = lO^'^ GeV. 

6 Conclusions 

We have investigated baryogenesis by R-parity violating squark decays in the frame- 
work of split supersymmetry. We have restricted ourselves to the baryon number 
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violating couplings Xijk to avoid rapid neutralino decay. These couplings involve 
only right-handed squarks. We assume that baryogenesis is dominated by the light- 
est one, which we take to be the right-handed up squark. A CP-asymmetry in the 
squark decays arises from the interference of tree-level and 2-loop diagrams involv- 
ing the CP-phases of Xijk- The relevant Boltzmann equations include, in addition 
to the baryon number violating interactions, also the supergauge interactions, the 
strong sphalerons and the top Yukawa interaction. The generated baryon asymme- 
try can be enhanced by some hierarchy in the Xijk, where the CP-asymmetry can 
be increased while keeping the washout small. 

The observed baryon asymmetry can be successfully generated if the up squark 
mass is larger than 10^^ GeV. A stronger constraint is induced by neutralino de- 
cays. In order to keep the lightest neutralino sufficiently stable to provide the dark 
matter, the up squark mass has to be at least 10^^ GeV. Depending on the gluino 
annihilation after the QCD phase transition, such large squark masses may induce a 
too large gluino lifetime. However, somewhat hghter left-handed squarks can speed 
up the gluino decay. At the same time the neutralino remains sufficiently stable 
and baryogenesis is hardly affected. Favorably, the LSP is bino-like with some hig- 
gsino admixture and a mass not far above while the gluino mass is in the 

few hundred GeV range. Because of the high sfermion mass scale, this scenario 
predicts sizable corrections to chargino and neutralino Yukawa couplings, which can 
be probed at a future linear collider [24] . 

It would be interesting to study also the case of the Icpton number violating 
couplings. These operators could be responsible for the light neutrino masses and one 
might wonder if this is compatible with successful baryogenesis. However, something 
would have to be added to the setup to replace the neutralino as the dark matter 
particle. 
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Appendix 



Using the approximations we discussed in sec. 3, the Boltzmann equations that 
govern the evolution of H, Yy, and Y^^ take the form 



dz 



(3K <j -Y_ + A{z 
+ [1-2B 



-^12 + r; , -^23 



5k + 2' 



'12 



12k + Q 
5k + 2 



B23)Yy + {B2s-B,2)Yd, 



d 
dz 



2k d 



-Y_ 



ik + 2 d 



-Y, 



5k + 2dz 5k + 2dz' 
= (3K['-^{Y^-Y^) + B,Y^-A{z)B,Yy} 



ma 352 



AioqCK.,— 1 



mi \ ^8k + 4: 



H Stt^'''''"^' z' 




m% j \5k + 2 



12fc + 6 

5k + 2 "^'^ 



^^]as-, I A?,, + ^?3AL3^ 1 (Y. + 2Yy + Y,,) 



H 37r2 z^ 



4 

m - 

as , 



^ 5k + 2dz 5k + 2dz ^ 



O- O 2^ 4 1 112 I "13"123 4 

•"''nA5Aif2y,-||±|y.-?^yJ. (44) 



The scatterings arc given for temperatures considerably smaller than m^, i.e. z 3> 1. 

In the special case B12 = -B23, the down number densities always appear in the 
combination Y^^ + TY^^ + l^s, since 7512 ,Ti2 — 1S23,T2z- (This is not true for the 
scatterings induced by o?3, which however are suppressed by m\lmS ~ 1/100 and 
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can be neglected.) Thus we can form a closed set of equations for YL and Yb 



dY_ 
dz 



(3K \ -y_ + A{z) 



Ik + 2 



Ik 



8k + 2 4k 



-B 



12 



5A; + 2 llA; 
+ 



%k + 2 4/c + l 



Bi2 I Yj 



B 



dY, 



B 



dz 



1^^) - 2BuY. 



5k + 2 Ilk + 5^^ 



4:k + l 



4:k + l 



- A{z)B,2 
171,-, 352 



2k 



4:k + l 
1 



6A'; 



4:k + l 



Yn 



H 37r' 



2-^112'^s^4- 



(45) 
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